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Abstract 

The motion of holes on the triangular lattice is studied using the t-J model. 
Within the Born self-consistent approximation and the exact Lanczos diag- 
onalization, the single hole physics is first analyzed. Then the spiral phase 
theory of Shraiman and Siggia is used to investigate the case of a finite density 
of holes. 
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I. INTRODUCTION 



The main motivation for the recent innumerable works on the Hubbard model is the 
Anderson suggestion^! that says that this model, in its simple version, can explain the physics 
of the CuO planes in high critical temperature superconductors. The Hamiltonian of this 
model is written as follows 

H = -tJ2( c l<T c j,<r + h.c \ + Uj2 n i,1 n i,l (1) 

where nj jCr = c\ a Ci )(7 is the occupation number of an electron with a spin a at the site i. 
The first term of H is the kinetic energy which allows for an electron to hop from one site 
to one of its nearest neighbors with an amplitude t. The second term stands for the on- 
site Coulomb repulsion (U > 0) and c| CT and c i)CT are respectively creation and annihilation 
operators. 

For large U, the double occupancy of a site with two electrons j and j is energetically 
discarded. In this limit the Hubbard model becomes equivalent to the so-called t-J model: 

H t -j = Po{-tE ( C W + h - c ) }^o + ^E ( S * • S i - l^rij) (2) 

(hi) y (m) 

where P = JI^ (1 — n^n^) applies the constraint of non double occupancy. S = 
\ YLafi c l,a a a,i3 c i,i3 * s the spin operator (a x , a v and o z are the Pauli matrices). The cou- 
pling constant of the magnetic part, J = 4t 2 /U, is positive. 

At half filling (one electron per site on average), the first term of H t -j is effectiveless 
and the model reduces to the Heisenberg model. It is now widely believed that the ground 
state (T = Oif) of the Heisenberg model on the square and triangular lattices presents 
antiferromagnetic long rang order in the thermodynamic limiti. 

In this paper, the motion of a single hole is first studied on the triangular lattice using the 
slave fermion representation and an exact diagonalization. Then the collective instabilities 
for a finite density of holes are investigated in the framework of a mean-field approach first 
introduced by Shraiman and Siggiai and known to lead to spiral phases. The main result 
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in that part of our work is that spiral phases in the doped case can only develop within the 
preexisting plane of the 120° AF structure of the pure system. The discussion of the present 
results is made together with a comparison with the square lattice results. 

In the second section, the slave fermion representation is used to write down an effective 
Hamiltonian from (fj). Section III describes the self-consistent Born approximation within 
which the single particle properties are calculated. The results, analyzed in sections IV 
and V, distinguish between the cases J ^ and J = 0. For the former case a Fermi 
liquid description can be done. However, for the later one and for positive t, the quasi- 
particle vanishes away from the centre of the Brillouin zone. The question of a 'liquid' 
where the quasi-particle depends on the position of the wavevector on the Brillouin zone 
arises naturally and will be addressed elsewhere. Section VI is devoted to an exact calculation 
using the Lanczoc method. Good agreement concerning the position of the energy minima 
as a function of k is obtained. Finaly, the motion of a finite density of holes is considered 
in sections VII and VIII where the spiral phase is calculated. 

II. THE MODEL 

A. Slave fermion representation 

The t-J model is characterized by the existence of spin and charge degrees of freedom. In 
one-dimension (ID), these degrees of freedom separate and one ends with a free propagation 
of the hole and the domain wall in the AF background. This illustrates the mechanism of 
spin-charge separation in ID. In 2D, a string of unsatisfied bonds takes place behind the 
moving hole. Spin and charge degrees remain coupled. 

The constrained electronic operator is written as follows 

c ii(T (l - n ijCT ) = c i;CT = iplb i>a (3) 

where ip\ is a fermionic operator which creates a spinless hole at the site i and b\ a is the 
Schwinger bosonS operator which creates a spin a at i. The number of bosons and fermions 
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on each site must satisfy the following constraint 

^ + E & Uv = 2S (4) 
The spin operator is written in the form 

Si = xl°Xi, (5) 

where 

is a two-component spinor. The Schwinger bosons formalism is easily generalized to higher 
values of the spin S where an expansion on 1/2S allows one to recover the spin wave theory, 
that is 



'2S - b\bi \ ( V2S - b\bi/2V2S \ 

, )"( * ) (7) 

for large S. Here 6j is the Holstein-Primakoff bosonic operator describing small deviations 
around the classical spin configuration which minimizes the energy (this configuration is 
hereafter called classical Neel state as on the square lattice): 

*=(7) (8) 

The constraint (P is transformed into a non holonomic constraint 

4^ + b\b t < 2S (9) 

which is more difficult to handle than (^). The latter is taken into account by introducing 
a site dependent Lagrange multiplier Aj and the effective action is written as follows 

dr <j -t ipl d T ipi + b\ a d T b if(y - t 



+ J2XiUHi + xlxi-2s'\\. (10) 

i 

In the spin wave theory the constraint is neglected and we do so in the present work. 



B. Effective Hamiltonian on the triangular lattice 



The classical Neel state on the triangular lattice is displayed in Fig. |T[ The spins are 
oriented, on a coplanar configuration, in such a way that two adjacent spins make an angle 
of 27r/3. These different orientations define three sublattices A, B and C instead of two 
sublattices as it is the case on the square lattice where two adjacent spins are antiparallel. 
By carrying out a local rotation of 27r/3 and — 2n/3 on the sublattices B and C respectively, 
that is 



Xb(c) = K B {c) ( j (11) 

where 1Zb(c) is the SU(2) matrix representing the rotation on B (C), the system is made 
formally ferromagnetic but the resulting Hamiltonian keeps naturally its AF character: 

R t _ J ^R=J- E(V>l% +hc) - f - + he) 

-^(S>te(fc - h \) - X>,Vi(fc - &})) 

+ V E 6 & + i E + - + • (12) 

The hole moves as a free particle (coherent motion) through the first term in (O) which 
is present since the spinors on adjacent sites are not orthogonal. In the case of the square 
lattice, such a term is absent. The second term in ( |12D allows an exchange of a spin flip 
and the hole by conserving the total number of overturned spins. In the third term, there 
is a change in the sign depending on wether the spin flip is absorbed or created after the 
hopping of the hole. J2' refers to the propagation of the hole in the three directions 
(— e$), Fig. |]. The last term is the AF exchange interaction written in spin wave theory. 

The presence of coherent and incoherent propagations of the hole is a salient feature 
of the triangular antiferromagnetic lattice, which in some sense interpolates between the 
ferromagnetic and antiferromagnetic limits on the square lattice. In the last case, only 
incoherent propagation is possible in the absence of quantum spin fluctuations. Arguments 
taken from the Brinkman-Rice picturei shows that the most important hopping terms in 



(|12|) are the first and the third. The point is that the coherent propagation alone would give 
a minimum in the hole dispersion of St for t > or of y for t < 0. But the hole can still 
lower its kinetic energy in a quite significant manner by using the incoherent channel. Indeed, 
would the third term be alone, then the estimate — t^/lE obtained in the self-retracing path 
approximation would apply. On the other hand, we expect the second term to play only a 
minor role and shall therefore neglect it in the following. 
Transcribed in k space, Eq. ( |i~2"D yields 

k,q 

~t 7k^k^k + JJ2 W q Q i«q ( 13 ) 



where 



and 



7 k = ^cos(k • e;), /i k = ^]sin(k ■ e,) 



The operator a q obtained from 6 q by using the Bogolubov transformation, creates a spin 
wave with a wavevector q. u q and t> q are the coherence factors: 



2cj q V 2 4 J 



t/ q = sign(-7 q j 



J 3 7„ LO, 



_ q 

2cj q V 2^4 J 



In the following, (|i~3|) is investigated using the self consistent Born approximation as devel- 
opped ini. 

III. SELF-CONSISTENT BORN APPROXIMATION 

The single particle Green's function is defined as follows 

G a (k,iu) = (Solck—L-Ck^o) (14) 

UJ — ±1 
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where |$o) is the ground state of the AF Heisenberg part of the Hamiltonian. Using ([|) one 

gets 

Cfc ' CT = "TtF 1 l>h b k+k',a (15) 

VN k , 

for the electronic operator in the slave fermion representation. As an approximation, the 
operator bk+k',a is replaced by its mean value in the Neel configuration since we consider that 
the important physics is contained in the fermionic part ip. We define the Green's function 
related to if) by 

G(k,u) = ($ \ip k —!—4\® ). (16) 
uj — tl 

In the Neel state, we have 

{bi,i) ~ exp (2im7r/3), and (b^i) ~ exp (—2imir/3) 

where m takes 0, 1 and —1 on the three sublattices A, B and C respectively. This can be 
written in the form 

(fe i>T ) « exp (— iG ■ Rj), and (b it i) m exp (iG ■ Rj) 
where G = (47r/3)x. Transcribed in Fourier space, these yield 

<6*,T> ~ V^4,-g, and (b k>l ) « VN6 k>G . 

A simple correspondance between the true Green's function, ([14]), and the ^-Green's func- 
tion, (|16|), namely : 

G T (k,cu) «G(-k-G,w) 

G ; (k,w) « G(-k + G,cj) (17) 

allows the evaluation of the true electronic Green's function G a using G. In this section, we 
concentrate on the Green's function of the slave fermion. 
The Dyson equation for G(\c,u) is 



where e(k) = —tj(l&) is the free part of kinetic energy and £(k, uS) is the self energy resulting 
from the incoherent motion of the hole. The evaluation of S(k, a;) is done using the diagram 
of Fig. |2| which corresponds to the following expression 

£(k, u) = —Y, |M(k, q) | 2 GW(k + q, u - W (q)) (19) 

q 

where is the ^-Green's function of the free part of the Hamiltonian. The self-consistent 
approximation consists in replacing in (|19"D by the total Green's function G. In this 
approximation, fllTf ) becomes 



/ It 2 \ - l 

G(k, u)=[u J - e(k) -— E \ M (^ q) l^(k + q, u> - w(q)) J (20) 

where M(k, q) = u ? /i fe - u q h k+q . 

In this approach, a series of an infinite number of diagrams which neglect vertex correc- 
tions is considered as shown in Fig. ^. Vertex corrections roughly correspond to Trugman 
processes^, which are known to become important fot t ^> J. In the most elementary process 
of this type on the triangular lattice, the hole makes one turn and a half around a triangular 
plaquette, which gives rise to coherent nearest-neighbor hopping. Trugman processes tend 
to push the energy minima towards the corners of the hexagonal Brillouin zone (independent 
of the sign of t). Their influence should remain week for moderate values of j-, as on the 
square lattice. 



IV. RESULTS 

Eq. (BO) is computed by means of an iterative method. We fixe an arbitrarily initial 
function G(k, u) for all wavevectors k, belonging to the hexagonal Brillouin zone, and fre- 
quencies lu, and use Q20D to iterate until convergence towards the unique solution. The 
unicity has been proved by starting from different initial functions. 
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The numerical calculation is performed for hexagonal clusters where the sites are labelled 
as shown in Fig. |^ and parametrized as follows : 

2tt 

k x = —m, < m < 3£ 
3£ 

4-7T / 1 1 - (-l) m \ 

to get the right number of independent sites inside the first Brillouin zone. The results 
presented here are from a cluster of 108 sites (I — 6). The number of independent sites is 
3£ 2 . 

The quantity of interest is the spectral function which is related to the Green's function 
through the relation 

A(k,u) = --lmG(k,u). (22) 

7T 

At finite J, -4(k, uj) presents a k-dependent quasi-particle peak at low energies. This is 
consistent with a Fermi liquid picture for positive and negative t. Note that because the 
electron-hole symmetry is absent for the t-J model on the triangular lattice, the spectral 
function depends on the sign of t. This is not the case on the square lattice where there is 
a t — > — t symmetry. 

For negative t and finite J, the energy minimum of the quasi-particle peak is located 
at k = (n,n/\^3) which corresponds to the middle of the edges of the Brillouin zone. 
For positive t, however the minimum is realized at the centre of the Brillouin zone: k = 
(0, 0). The spectral functions of these wavevectors, together with those of k = (4tt/3, 0), are 
displayed in Fig. |] and Fig. |] for J/\t\ = 0.2. 

V. PROPERTIES OF THE QUASI-PARTICLE 

A. Spectral weight of the quasi-particle for t < 

In this section, we are interested in analyzing the properties of the quasi-particle in the 
case of negative t. The spectral weight a(k) of the quasiparticle is calculated by computing 
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the area under the quasiparticle peak. The results for small values of J suggest a simple 
power law as a function of J: 

a(k) ~ J a (23) 

where the values of a and (3 are summarized in table I for the three characteristic wavevectors 
k= (0,0), (4vr/3,0) and 

Contrary to the case of free electrons where the weight is a constant (= 1), a(k) is here k- 
dependent. The zero quasi-particle weight at J = implies the absence of the quasiparticle 
peak at low energies as shown by the spectral function. What the absence of quasi-particle 
means is the breakdown of the Fermi liquid picture. However, Nagaoka theorem! applies for 
negative t and ensures that the ground state is ferromagnetic at U = oo (J = At 2 /U). Our 
result in this case is an artifact of the approximation which states an AF background. For 
finite J, it is natural to consider such a background and our results are physically meaningful. 

For J ^> |t|, the spectral weight goes to unity since 

1 



a ( k ) = gTFTTwTI ~ , oo = E(k) 



l-£ReE(k,u;) 



and 



|M(k ' q)P 



Olo v ' ' ' N ^(^(q)-t 7 k) 2 J 2 

B. quasi-particle dispersion relation for t < 

By performing the numerical calculation of the position of the quasi-particle peak for 
every wavevector belonging to the Brillouin zone, we get the dispersion in k of the energy 
minimum. This energy is well fitted by the following simple expression 

E(k) =A + B lk + CK 2 (24) 

suggested by the arguments developped in!. A, B and C are J-dependent as seen in table 
II. 
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This fit contains hopping processes to nearest neighbors through the 5-term and second 
nearest neighbor through the C-term. comes from the coherent part of the Hamiltonian, 
whereas Ch^ 2 originates from the incoherent part with the partcipation of quantum spin 
fluctuation. So at finite small J, a Fermi liquid picture is appropriate to describe the motion 
of the hole as a quasi-particle with the dispersion relation given by 

For J ^> \t\, one gets an analytical expression for E(k), namely 

B(k) ._ i7k _ £ J«l. (25) 



C. The case of t > 

At finite J and positive t, the physics is similar to that of the case of t < 0. But for 
J = 0, the situation is quite different since the Nagaoka theorem is not satisfied and the 
energy of the hole can be minimized further in a singlet spin state. It is then natural to 
consider the AF background as a good approximation for t > even when J=0. From the 
calculation of the spectral function .A(k), we found out that the energy minimum is located 
at k = (0, 0). For k = (0, 0) ( J = 0), a well defined quasi-particle peak is present. However, 
the quasi-particle is strongly k-dependent. the peak loses in intensity and broadnes, as 
illustrated in Fig. [5], when we run away from the center of the Brillouin zone. Therefore, 
depending on the wavevector, a Fermi liquid interpretation, in the classical way, is either 
possible or not. We believe that Trugman processes will not spoil this conclusion since they 
only renormalize (slightly) the coherent part of the hole motion. This discussion rises the 
question of an electronic system where the quasi-particle weight depends on the wavevector 
k and vanishes on some points of the Brillouin zone0. A work concerning this question is 
in progress. 



11 



VI. EXACT DIAGONALIZATION 



The exact diagonalization is performed using the Lanczos method on a 3 x 2 2 hexagonal 
cluster (Fig. |3|). The dispersion relation and the total spin are calculated. 

The results for t < and t > are presented. For t < the dispersion relation E(k) 
and the total spin S to t are summarized in tables III and IV respectively. E(k) is reported 
as a function of the wavevector k and J/\t\. First let us consider the case t < 0. For 
< J < 0.3, the energies have a small dependence on k. But for J > 0.3, the energy 
minimum is located at k = (27r/3, 0) which transforms into the middle of the Brillouin zone 
sides, namely k = (tt, 7r/v3), by a translation G = (47r/3,0). This is consistent with what 
is obtained in the approach of the previous sections. However, a translation of — G produces 
a meaningless result. The reason for this discrepancy is due to the fact that in the slave 
fermion approach the chiral symmetry is broken, whereas in the exact diagonalization the 
symmetry cannot be spontaneously broken. 

Things become more understandable if we calculate the total spin S to t- What is clear 
from the table IV is that for J > 0.3, Stot is small for all the values of k and the ground 
state is a singlet: S tot = 1/2. For < J < 0.3, Stat is big: 5/2 < S tot < 11/2. For J = 0, the 
ground state, given at k = (0, 0), is ordered ferromagnetically, a result which agrees with the 
Nagaoka theorem which applies only to negative t on the triangular lattice. For < J < 0.3 
and k ^ (0,0), the energies are very close to —6 (the lowest energy). The ferromagnetic 
states have, however the energies 2, —1 and 3 (in the units of \t\) for (71,11/^), (2n/3, 0) 
and (47r/3, 0) respectively. A possible explanation is that the system finds a compromise 
in which the spins deviate slightly from the ferromagnetic state (the spins do not remain 
in coplanar positions) to keep a high value of Stat but the energies are close to —6. This 
phenomenon is similar to Aharonov-Bohm effectS. The stability of this phase as the size 
of the cluster increases has to be clarified. On the other hand, it is natural to think that 
the AF correlations, growing with the cluster size, will rise above any other correlations at 
J > 0. 
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For large J, the antiferromagnetic correlations become dominant and the total spin 
becomes small. In this case the exact results of the minimum of the dispersion relation 
compare well with the self-consistent approach. 

The situation is simpler for positive t. The results are summarized in tables V and VI. 
There is no crossing in the energy levels as it is the case for t < (table III). The minimum 
is always located at k = (47r/3, 0) which transforms into k = (0,0) using Eq. fll7|) . The 
total spin of the ground state is S to t = 1/2. Here also a good agreement, concerning the 
minimums of EQ&), between the two approaches is obtained. For J = 0, the fact that the 
ground state is not ferromagnetic confirms that the Nagaoka theorem does not apply for 
t > 0. 



VII. MOTION OF A FINITE DENSITY OF HOLES 



We assume that the holes live around the minima calculated in the case of one hole. 
They form small pockets or valleys, whose area is proportional to their density. We shall 
use the slave fermion picture from now on. There is only one valley around the centre of 
the Brillouin zone for t > but three of them on the edges of the Brillouin zone for t < 0. 
We are interested in the long-range interaction between holes, as mediated by low energy 
spin- waves. 

The expression of the coupling of holes to the spin waves is 



- H— ^I+q^k (Uq/lk - ^q/ik+q)aq + (V&k ~ Mq^k- 



k,q 



.qjQ!_q 



(26) 



Among the three Goldstone modes (q = and q = ±(47r/3, 0)), only q = is relevant since 
it is easily seen that two valleys cannot be coupled by a momentum transfer q = ±(47r/3, 0), 
Fig. |] in the case t < 0. This is a fortiori true in the case t > where there is only one 
valley left. For q ~ 0, ( PS] ) becomes 

\/3 r 



k,q 
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-^T^S^k+qV'k {u q + v q ){-q.Vh k ){a q + al q ) (27) 



AN . 

k,q 



which transforms into 



- ^E4qtt(-^ z (g)) + ^E4 + qV>k(p*.(zq0(q))) (28) 

k,q k,q 

where 

S*(q) =<(&„- &L,)/2 (29) 

is the slowly varying component of the magnetization in the z-direction, while 

0(q) = (6 9 + 6 t _ g )/2 (30) 

parametrizes the slow distortion of the ordered 120° sructure within its plane. The vectorial 
quantity 

Pfc = Vfc/ifc (31) 

appears like the effective dipolar momentum carried by the hole. It can be shown that 
the coupling to S z (q), which semiclassically behaves as does not lead to long-range 
interaction between holes. On the other hand, the second term, once written in real space, 
becomes 

^EE^WV,(r)(p,V r 0(r)) (32) 

r i 

where the i indicates the different valleys (see Fig. |6|) . The effective dipolar momentum can 
safely be taken equal to its value at the centre of each valley for small hole density. Two 
cases are to be considered: 

i) t > 0: the minimum of the dispersion relation is located at k = 0. There is only one 
valley and the dipolar momentum p = 0. So, no dipolar interaction between the holes and 
the spin waves can be generated. The interaction is at least quadrupolar and decays as r~ 4 . 
it) t < 0: the minima of the dispersion relation are located at k = (tt, and equivalent 
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momenta, Fig. |6|. The dipolar momenta in this case do not vanish. They take the following 
values 

Pi = 2x, p 2 = -x + v^y, p 3 = -x - v^y. (33) 

VIII. SPIRAL PHASES 

According to the results of the previous sections, we introduce a phenomenological Hamil- 
tonian for the motion of a finite density of holes. The spatial density of this Hamiltonian is 
given by: 

f 1 11 M 2 

n = E " ^JV 2 ^ + <?aV>M(pV)0 + -(p(V0) 2 + — ) (34) 

j=l,2,3 L A J ^ X 

where M is the z-component of the magnetization (M ~ 0). The coupling constant g is of the 
order of t for t <C J {g = tay/3/4). When J <^ t, the vertex corrections become important 
and are expected to renormalize g to an order of J. Using the expressions of Pi,P2,P3 
given above, it can be seen that d x <j) and d y <j) are respectively coupled to 2n\ — n 2 — n 3 and 
V3(n 2 — Tis) where rii = (ipjipi) is the hole density in the valley i. At the mean-field level, 
we can write: 

n in t = sa(Re{m + n 2 e~ 2i ^ 3 + n 3 e 2 ^ 3 )(«9 x .0 + id v <j>)\ + \d x 4> + id y( j)\ 2 '. (35) 
The minimization with respect to <p yields 

d x <f> + id y <f> = -j(ra + n 2 e~ 2 ^ 3 + n 3 e 2 ^ /3 ). (36) 

Eq. ( |3~5"1) becomes 

1 n 2 n 2 

1 y a „-2mt/3 , _ „2i7r/3|2 /o 7 \ 

riint — ~- \n 1 + n 2 e + n 3 e \ . {6() 

Z p 

and the kinetic energy can be written as follows 

+ n l + n l ) = JL ({n x + 7l2 + n ^ 2 + 2\n 1 + n 2 e~ 2 ^ 3 + n 3 e 2i ^ 3 \ 2 ) (38) 
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so that the normal phase becomes unstable when 

2tt 1 g 2 a 2 



< 



3m 2 p 

that is when 3g 2 a 2 m/ pn. Once this condition is satisfied, the system maximizes 
|m + n 2 e- 2i7T/3 + n 3 e 2i7r/3 \ 2 = I ({n x - n 2 ) 2 + (n 2 - n 3 ) 2 + (n 3 - n,) 2 



2 

This is realized for n\ = n > and n 2 = ^3 = or all the solution obtained by cyclic 
permutation of the indices 1, 2 and 3. This implies that only one valley is occupied and the 
two others are empty in the spiral phase. A simple interpretation of this phase in the real 
space is obtained using Eq. (j3~6|). For the phase ri\ = n and n 2 = n 3 = one has 

(d x <j)) = -^ n , an d (d y <f>) = 0, 

which means that the spins rotate around their position in the normal phase uniformly when 
moving along the x-axis. In the other solutions for rij, the rotation occurs along the axes e 2 
and e 3 for n\ = n 3 = and n 2 = n, and n\ = n 2 = and n 3 = n respectively 
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FIGURES 

FIG. 1. The classical Neel state in the case of the triangular lattice is drawn on one plaquette. 
The three sublattices A, B and C are shown as well, correspond respectively to the vectors on 
A, B and C sublattices for i = 1, i = 2 and i = 3. 

FIG. 2. On the left hand side of the figure, the diagram used for self-energy is shown (with the 
bare propagator substituted by the true one). The expansion in the bare propagator is shown on 
the right hand side. The wavy line refers to the spin-wave propagator. 

FIG. 3. The Brillouin zone is displayed for 1 = 2. In this case, the wavevectors k = (0,0), 
(n,ir /3 l l 2 ) and (47r/3, 0) correspond respectively to the points (m,n) = (0,0), (3,0) and (4,0). 

FIG. 4. The spectral functions of: (a) (0,0), (b) (47r/3,0) and (c) (ir,ir /3 1 / 2 ) are presented as 
a function of frequency. I = 6, t = — 1 and J = 0.2|t|. 

FIG. 5. The spectral functions of: (a) (0,0), (b) (47r/9, 0) and (c) (7r,7r/3 1//2 ) are presented as 
a function of frequency. £ = 6, t = — 1 and J = 0. k = (4-7r/9, 0) is located near the point (2,0) of 
Fig. 3. The quasi-particle broadens. 

FIG. 6. The position of the valleys and their dipolar momenta pi are presented on the Brillouin 
zone for i = 1, i = 2 and i = 3. 
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TABLES 

TABLE I. The values of the coefficients a, and (3 of Eq. (23). 
TABLE II. The values of the coefficients A, B and C of Eq. (24). 
TABLE III. The hole energy as a function of k and J for negative t(= — 1) 

TABLE IV. Total spin as a function of k and J for negative t(= —1). 
TABLE V. The hole energy as a function of k and J for positive t(= +1). 
TABLE VI. Total spin as a function of k and J for positive t(= +1). 
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k 


P 


a 


(0,0) 


0.112 


1.056 


(4tt/3,0) 


0.440 


1.473 


(7T,7T/V3) 


0.517 


0.610 



Table I 
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J 


A 


B 


C 


0.05 


-4.186 


0.019 


0.013 


0.1 


-4.018 


0.032 


0.025 


0.2 


-3.766 


0.049 


0.040 


0.3 


-3.585 


0.070 


0.058 


0.5 


-3.271 


0.107 


0.077 



Table II 
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J k 


(0,0) 


(n,n/V3) 


(2tt/3,0) 


(4tt/3,0) 


0.00 


-6.000 


-5.410 


-5.662 


-5.348 


0.05 


-5.258 


-4.856 


-5.040 


-4.807 


0.10 


-4.517 


-4.302 


-4.418 


-4.276 


0.20 


-3.368 


-3.342 


-3.349 


-3.382 


0.30 


-2.917 


-2.907 


-2.898 


-2.838 


0.50 


-2.368 


-2.195 


-2.454 


-2.065 


0.70 


-1.832 


-1.562 


-2.073 


-1.609 


0.80 


-1.567 


-1.359 


-1.895 


-1.483 



Table III 
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Jk 


(0,0) 


(7T,7T/VS) 


(2tt/3,0) 


(4tt/3,0) 


0.00 


11/2 


9/2 


9/2 


7/2 


0.05 


11/2 


9/2 


9/2 


7/2 


0.10 


11/2 


9/2 


9/2 


7/2 


0.20 


7/2 


5/2 


5/2 


5/2 


0.30 


1/2 


3/2 


1/2 


3/2 


0.50 


1/2 


3/2 


1/2 


3/2 



Table IV 
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J k 


(0,0) 


(7T,7r/>/3) 


(2tt/3,0) 


(4tt/3,0) 


0.00 


-4.230 


-4.120 


-4.149 


-4.270 


0.05 


-4.076 


-3.928 


-3.971 


-4.192 


0.10 


-3.941 


-3.760 


-3.810 


-4.117 


0.50 


-2.964 


-2.706 


-2.699 


-3.569 


2.00 


0.574 


0.208 


0.328 


-1.725 


15.0 


18.068 


16.615 


16.383 


13.496 



Table V 
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Jk 


(0,0) 


(7T,7T/VS) 


(2tt/3,0) 


(4tt/3,0) 


0.00 


3/2 


1/2 


3/2 


1/2 


0.10 


3/2 


1/2 


1/2 


1/2 



Table VI 
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